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Abstract. Theoretical study of the indirect coupling of nuclear spins (qubits) 
embedded into a mesoscopic ring and in a finite length quantum wire in a magnetic field 
is presented. It is found that the hyperfine interaction, via the conduction electrons, 
between nuclear spins exhibits sharp maxima as function of the magnetic field and 
nuclear spin positions. This phenomenon can be used for manipulation of qubits with 
almost atomic precision. Experimental feasibility and implications for quantum logics 
devices is discussed. 
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Introduction 

The possibility that a computer could be built employing the laws of quantum 
mechanics has stimulated huge interest in searching for useful algoritms and realizable 
physical implementations pp. There are currently many promising approaches to 
quantum computation the most promising solid state approaches are based on the 
superconductivity 0, conduction electron spins |1] and nuclear spins 001111313110] as 
qubits. The attraction of the approach of using nuclear spin in quantum computation 
lies in the idea of incorporating nuclear spins into a semiconductor device EI. 
Using of the nuclear spins incorporated in a heterostructure in the highly nondissipative 
quantum Hall effect regime, as qubits [Sj, are promising because they are extremely well 
isolated from their environment and have a long decoherence time ^] . The 

energy gap in the spectrum of two - dimensional electrons in a strong magnetic field 
imposes severe restriction on the flip-flop processes, since the electron Zeeman energy 
is orders of magnitude larger than the nuclear one It follows that the nuclear 
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spin imbedded in 2DES under the QHE conditions is practically decoupled from the 
conduction electron spins. As a result the nuclear spin relaxation time has an activation 
behavior, i.e. is exponential in the electron energy gap and inverse temperature ^Hj 
and can be manipulated by external magnetic field and sample parameters over several 
orders of magnitude time interval. 

A set of five essential criteria for the physical realization of a quantum computer 
was formulated by DiVincenzo and coworkers [I4j. They are: 

1. A scalable physical system with well characterized qubits. The system should 
consist of a collection of independent subsystems each with a two-dimensional Hilbert 
space, so-called quantum bits, or qubits. Its physical parameters should be accurately 
known, including the internal Hamiltonian of the qubit, the presence of and couplings 
to other states of the qubit, the interactions with other qubits, and the couplings to 
external fields that might be used to manipulate the state of the qubit. 

2. The ability to initialize the state of the qubits to an initial state. This arises 
from computing requirement that registers should be initialized to a known value before 
the start of computation. Moreover, the initialization of qubits is used in quantum error 
correction algorithms. 

3. Long relevant decoherence times, much longer than the gate operation time. It 
was shown that for fault-tolerant quantum computation the magnitude of decoherence 
time scales should be 10^ — 10^ times the "clock time" of the quantum computer, that 
is, the time for the execution of an individual quantum gate. 

4. A 'universal' set of quantum gates, a) It should be possible to perform precise 
unitary operations on the individual qubits. b) Furthermore, the inter-qubit interaction 
should be controlled with almost atomic precision. 

5. A qubit-specific measurements capability. Any quantum computer would deliver 
its output as results of measurements performed on it. 

In what follows we concentrate on nuclear spin based solid state models El 13 
[TT| ISl ^1 inj. All the existing models satisfy the criteria 1,2 and 4a. In some models 
the criteria 3 is almost satisfied. At least, the ratio of gate to decoherence time allows 
to create a few qubit quantum computer. No existing model addresses the problem of 
controlling the inter-qubit interaction with atomic precision. In this paper we propose 
a new system in which criteria 1-4 are completely fulfilled. Our investigation is mainly 
focused on the criterion 4b. 

The proposed system consists of nuclear spins (qubits) embedded into a zero nuclear 
spin mesoscopic ring or a finite length quantum wire (Fig. la and Fig. 2a). The 
hyperfine interaction of the electrons in the system with nuclear spins leads to an 
effective indirect nuclear spin interaction. In what follows we calculate the effective 
nuclear spin interaction energy. The effective nuclear spin interaction energy can be 
chosen to have the following form: 




(1) 



where /j is magnetic moment of a nucleus, A and B are functions of the system 
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parameters as described below. We obtain that the effective nuclear spins interaction 
exhibits sharp maxima as function of the magnetic field and nuclear spin positions which 
opens the way to manipulate qubits with almost atomic precision. The selective nuclear 
spin interaction can be obtained by changing external parameters of the system. 

The first calculation of the indirect, via electron spins, hyperfine coupling between 
nuclear spins, was performed by Ruderman and Kittel J3] for the case of 3D metal in 
the absence of magnetic field (see pTOJ). Influence of mesoscopic effects on the RKKY 
interaction was studied by Spivak and Zyuzin ^^l- They pointed out that the only 
difference between pure and impure metals is an additional random phase which depends 
on impurity distribution. A mechanism of the indirect, via the exchange of virtual 
electron-hole pairs (spin excitons), nuclear spin interaction in the quantum Hall effect 
systems was suggested in [iHI and further elaborated in [Hj. Quantum computation and 
communication devices, based on this mechanism are proposed in [3 El IBjj 

Theoretical framework and results 

Let us consider a system consisting of electrons confined by a potential V (r) and 
interacting with two nuclear spins. We assume that the nuclear spins are located far 
enough from each other, so that the direct (dipole-dipole) nuclear spin interaction is 
negligibly small as is in isotopically engineered Si/Ge hetero junctions p. The contact 
hyperfine interaction between electrons and nuclear spins leads to an indirect nuclear 
spin interaction. 

The wave function of the electron 0(r) can be written as a product of an envelope 
function \l/(r) by the rapidly varying periodic function Mo(r) [1^ 

0(r) = ^{r)u',{r) = ^{r)uo{r) {V/ny^ (2) 

where uo{r) is the k = Bloch state, V is the volume of the sample, and fl is the volume 
of the unit cell. The function UQ{r) is conveniently normalized in the cell volume Q, 

J \uQ{r)f dr = 1. The localization of electrons is described by the envelope part of the 

n 

wave function \E'(r). The norm of the envelope function is J {"if {r)f dr = fl [TH] . 

V 

In what follows we will consider the system in the envelope function approximation, 
the rapidly varying function Mo(r) appears only in the expression for hyperfine 
interaction. The Hamiltonian of the system is given by H = Ho + h[^^ + h[^^ with 

Ho = ^ (p+^a)' + V{r)- gf^saH (3) 

and 

H? = y/iB^7n |wo(r.)|' ahS (r - r^) (4) 

where Hq is the Hamiltonian of the electron in the mesoscopic structure in the magnetic 
field, Hi'^ + h[^^ is the perturbation due to the contact hyperfine interaction, m* is the 
effective electron mass, A is the vector-potential, g is the electron (yf-factor, /i^ is the 



Nuclear-spin qubits interaction in mesoscopic wires and rings 



4 



Bohr magneton, H is the magnetic field, 7„ is the nuclear gyromagnetic ratio, Ij and 
a are nuclear and electron spins, Fj is radius vector of i nucleus, i = 1,2. Because the 
electron wave function is strongly peaked on the nuclei, the contact hyperfine interaction 
energy greatly exceeds dipolar spin interactions. 

The effective nuclear spin interaction energy calculated in a second-order 
perturbation method is given by the expression jl6j : 



r(2) 



E=T. ^ J' 'J ^/^ (1 - //) + CO. (5) 



E,,Ef Ef - 

Here and Ei are eigenf unctions and eigenvalues of Hq and fi is the electron 
distribution function. In this Letter we restrict ourselves to the single electron 
approximation, which proved to be sufficient for clean low density quantum wires 
and rings, see [20] and references therein. Presence of impurities violates the Kohn 
theorem [21] and the electron interactions may play an important role. In strongly 
correlated dense ID electron systems the Fermi hquid approach should be replaced by 
the Tomonaga-Luttinger theory, see |22] and references therein. The infiuence of the 
electron correlations on the results obtained here will be the subject of a more detailed 
publication. 

Nuclear-spin interaction in mesoscopic rings 

Consider a torus-shaped quantum ring of inner radius a, thickness d <^ a and 
negligible height h |2S] in a uniform parallel to the axis of the ring magnetic field 
H in z direction with two nuclear spins located at r(p, z) = ri (a + d/2, ipi, h/2) and 
r = r2 (a + d/2, ip2i h/2) (Fig. la). In this subsection we use the polar coordinates. The 
electron confining potential V (p, z) is 

{0 ii a < p < a + d and Q < z < h , 
oo otherwise. 

Firstly, let us find eigenfunctions and eigenvalues of the Shrodinger equation with 
the Hamiltonian Q. Due to the axial symmetry of the ring, the wave function can be 
written as follows 




^™,n,.=± = \h^\\ , , , e^"^^i?™,„ (p) Z,{z) (7) 



We assume that in z direction the electron is always on the ground level of the one- 
dimensional quantum well of thickness /i, thus Zi[z) = ^J^sin (x)- O^^^ region 
a < p < a + d the radial part of the wave function Rm,n (p) satisfies the equation: 



2m* pdp dp ' 2mep^ \ $c 
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where $o is the magnetic flux quantum and the radial number n = 1,2,3,... . The 
boundary conditions imposed by the potential V (p, z) (0) are: Rm,n (o,) = and 
Rm,n (a + d) = 0. 

We will look for the solution of Eq.(jHJ assuming that aX d <^ a the vector-potential 
varies slow in the ring, so we can put $(p) ^ $(a + (i/2). In this case the solution of 
Eq. (jHJ is written in terms of Bessel functions: 

Rm,n (p) = CiJa {pam,n) + C2YA {pam,n) (9) 

and am,n = y — "^a™'" , the constants Ci, C2, energy 
levels Em,n are defined by boundary conditions and the normalization condition 

a+d 

I P^mn (p) -Rm,n (p) dp = 1. Unfortunately, the wave function of the form (jH)) allows 

a ' 

to calculate the effective nuclear spin interaction only numerically. 

To obtain an analytical result, consider the effective interaction energy constants 
A and B for the case of infinitely narrow ring. To do this, let us set p = a in (jH)). It is 
readily seen that in this case the radial part of wave function decouples from the orbital 
part, Rn{p) = \l\sin ^ (''-^)^" ^ and 

/7m\^ ( ^a + d/2)Y , , , 

= 2^;^ (t) + W + j ^ (1°) 

For the infinitely narrow ring it's reasonable to consider the states only with n = 1. 
With Eq.(0), the effective nuclear spin interaction constants are: 

A=K.Y. " 'r " - /(g...)) . (11) 

m,n -C'l,n,+ i^\,m — 

where Kr = 2 (^^^pBhin |no(0)|')'. 

Using the wave function in the form Eq.((Z)), we numerically calculate the magnetic 
field dependencies of the nuclear spin interaction constants A and B using the following 
set of parameters: y^i = 0, <y92 = ^r, a = lOOnm, h = Inm, d = Inm and 0.5nm, 
one electron in the ring, T = . The results of our calculations are shown on Fig. lb. 
The magnetic field dependence of the effective nuclear spin interaction constants are 
defined by the energy level statistics. The nuclear spin interaction constant A describes 
nuclear spin fiip-fiop processes which are performed through flips of electron spin, and, 
therefore, the main contribution to this process is due to the energy levels with the same 
orbital quantum numbers and different spin directions, what gives 1/H dependence of 
the interaction amplitude. The nuclear spin coupling constant B is connected with 
electron transitions between energy levels with different orbital quantum numbers, but 
the same electron spin direction. Their periodicity with the magnetic field results in the 
periodicity of the constant B. It is seen that at the values of the magnetic field, when 



where A 



m + 



<i>{a+d/2) 
<I>0 
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there are two ground states, function B has discontinuities. In this case the perturbation 
theory is not apphcable. The obtained result is in a good agreement with the result 
obtained for the case of infinitely narrow ring potential Eqs. ()llfl2|l . 



Nuclear-spin interaction in mesoscopic wires 

The next system under consideration is a finite length quantum wire of the length / in 
x-direction, of the thickness d and of the height h with two nuclear spins located at 
t(x, y, z) = Ti (Zi, d/2, h/2) and r = r2 {I2, d/2, h/2) in a magnetic field H in z direction 
(Fig. 2a). We suppose that the transversal sizes of quantum wire are much smaller 
than the length of the quantum wire and the cyclotron orbit of electron. We consider a 
model, when the confining potential is 

- - , -y- , - - ^^g^ 
otherwise. 

The eigenfunctions and eigenvalues of the Hamiltonian Q with potenital ()13|) are 




h'^n'^ 2 ^^^^ 2 ^^^^ ,2 / / \ 

En,m,k,^ = + + T 9t^BH/2, (15) 

As in the case of the nuclear spin interaction in mesoscopic wires, we assume that 
d,h <^ I, so we consider only the electrons on the ground levels of potential wells in y 
and z directions, i.e. m = k = 1. Substitution of Eqs. (fT^ and (fT^ into (0) gives us 
the effective nuclear spin interaction constants. 

We have calculated analytically and numerically the effective interaction constants 
A and i? as a functions of nuclear spin positions for odd and even number of electrons 
in the wire. The results of calculation for the odd number of electrons (A^ = 9) at T = 
are presented at Figures 2b-2d. In low magnetic field region, when Zeeman splitting 
energy is much less than the energy gap between levels with different n (Eq. ()15|) ) and 
at T = limit the expression for A takes a simple fortm: 

,, (^"'(^)-''(^))\ 



1 2 



2 



where K,„ = 2\ ^^^"^'^^2"^"°^°^^ 1 . In this limit the interaction constant A has a set of 



1 ■ildh 

^^^^^ maximums (Fig. 2b) and B <^ A. Fig. 2c shows an increasing of interaction 
constant A if nuclear spins are located not far from each other (in the vicinity of the line 
li = I2) and a decreasing of interaction constant A for the other nuclear spin positions 
with increasing of the magnetic field. The interaction constant B has a non-trivial 
dependence on the nuclear spin positions (Fig. 2d). 
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Conclusions and discussion 

To conclude, we proposed a new possible implementation of a basic unit for quantum 
computer based on the nuclear-spin qubits embedded into the zero nuclear spin 
mesoscopic ring or finite length quantum wire. Particular emphasis has been placed on 
the investigation of nuclear spin interaction via the electrons confined in such systems. It 
was found that the indirect nuclear-spin qubit interaction is very sensitive to the system 
parameters: nuclear spin location, number of electrons, magnetic field and geometry 
of the system. It dependence on the system parameters is completely different from 
indirect nuclear spin interaction in 2D and 3D metals. Preliminary finite temperature 
calculation indicates that the values of the effective nuclear spin interaction constants 
are decreasing with temperature, but the main features of the obtained results remain 
quahtatively unchanged. 

Now let us consider how our model of quantum computer is matched by the set of 
DiVincenzo criteria listed in the Introduction. 

1. It's well known that nuclear spins are appropriate candidates to be qubits piTUj. 
Nuclei with spin 1/2 are two-level systems with well defined states |0 > and |1 >. 
Needless to say that such a system is scalable and, basically, there is no principal 
limitation on the reasonable number of nuclear-spin qubits which can be integrated into 
a quantum circuit. The ring architecture of the quantum computer is of considerable 
promise. Let us imagine a mesoscopic ring with nuclear-spin qubits located in the 
immediate vicinity of it. A local change of the electrostatic potential near a qubit by 
a gate electrode changes the value of the envelope wave function on the qubit and, 
correspondingly, allows to switch on the interaction between any two qubits, whereas 
almost at all quantum computer proposal only adjacent qubits can directly interact. 
Nuclear-spin qubit interaction through electrons confined on a sphere opens a further 
way to improve the quantum computer architecture. 

2. We propose to initialize the nuclear-spin qubits using spin-polarized electrons. 
There is a general agreement that it is the only feasible method since the nuclear spins 
are highly isolated from the environment. Possibility of nuclear spin polarization by spin- 
polarized transport was demonstrated almost ten years ago [211 • Possible methods 
to introduce nonequilibrium polarizations of the electrons include injections of spins 
from ferromagnetic contacts [2^1 , optical pumping [2Z1I2H1 and spin refrigeration (see [7] 
and references therein). 

3. The decoherence time of the nuclear spins in mesoscopic systems is expected 
to be long enough to perform the quantum computation, since the discrete electron 
spectrum in mesoscopic systems imposes restriction on the fiip-fiop processes and the 
nuclear spin relaxation time at law temperatures is expected to have an activation 
behavior ^31- Up to the present, a calculation of nuclear spin decoherence time in 
mesoscopic structures is not made, however as a rough estimation we can take the 
decoherence time of the nuclear-spin qubit embedded into 2DEG, T2 = 10 sec jTHj. The 
characteristic time of qubit interaction is Tint ^ h/A, where for the nuclear-spin qubits 
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in finite length quantum wire A is given by Eq. Using available experimental data 

for GaAs, |mo(0)|75As ~ 9-^ ' lO^^cm"^ and following set of parameters: a = 200nm, 
d = h = 5nm, H = O.OIT we obtain Tint = 2-10~^sec and T2/Tint = 5-10^. Our estimate 
indicates that the present system is suitable for fault-tolerant quantum computation. 

4. a) The one-qubit operations using the NMR could be similar to the existing 
experimental suggestions |H1 E]. b) Interaction between any two qubits, which is 
necessary for two-qubit operations, is performed by the confined electrons. By varying 
the external parameters (magnetic field, number of electrons, gate potentials) we can 
control with almost atomic precision the nuclear spin interaction strength by creating 
maxima of the amplitude of electron wave function on some qubits and zero on the other. 
As an example, consider how it is possible to control nuclear-spin qubit interaction in 
the finite length wire by changing the number of electrons. Let us place the qubits 
at some predefined positions with coordinates Xi = where di are different prime 

Pi 

numbers and are integers (oj < f3i). The orbital part of the last electron level wave 
function (from Eq. ()14|l ) for 

M 

HA 

= r-f^ - 1, (17) 

electrons in the system is 



M ^ 

n A 

V J 



sm I — — I sin I 



It is readily seen that '^n,m,k,±{xi) = for i 7^ j, k. This means, that all qubits except 
j and k are at the nodes of the wave function and only j and k qubits interact. Fig. 3 
illustrates this idea. It follows that the accuracy of nuclear spin positioning should be 
few atomic units. 

5. In mesoscopic systems the single nuclear-spin measurement is still an open 
problem. One of the possibilities is to use Hyperfine Aharonov-Bohm effect (HABE) 
|29j . following from the coupling of the nuclear spin polarization to the phase of the 
conduction electron wave function. This was outlined in jH] and will be considered 
in details elsewhere. Another possibility is to use spin- dependent magneto-transport 
tunneling through a system, which energy levels are spin-split by an external magnetic 
field pni- The tunneling current develops a distinctive peak at the frequency of Zeeman 
splitting, which can be sufficiently narrow to measure a state of the nuclear spin. 
Moreover, it's possible to use spin to charge conversion as it was discussed in 
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Figure 1. a) Two nuclear spins embedded into a torus shaped quantum ring; 
b) dependence of nuclear spin coupling constants A and B on the magnetic field, 
P — 
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Figure 2. a) Two nuclear spins embedded into a finite length quantum wire in 
a magnetic field; b) nuclear spin coupling constant A in low magnetic field region 
j^igt^m^tjl _ Q iy^ nuclear spin coupling constant A in high magnetic field region 
(^iMs^l^LJL — 3) and d) nuclear spin coupling constant B on the nuclear spin positions. 
Odd number of electrons, T ~ 0, the nuclear spin coupling constants are given in units 
of Kyj/E^, where = ■ 
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Figure 3. Qubit arrangement and the last occupied electron level wave function for 
N = 3,5,9 electrons in the wire. The qubits are located in the nodes of the different 
wave functions. 



